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1 . Introduction 

Following the rapid growth of the helds of spintronics and molecular electronics [1, 2, 3], 
there have been great achievements in modern and condensed matter physics. These 
new and exciting helds of physics whose interests for are manifold, help in explaining the 
mechanism of spin hip transitions and spin propagations [4] in single- and multi- molecule 
magnets[5, 6, 7, 8], quantum computation [9, 10] and Landau-Zener(LZ)[ll, 12, 13, 14] 
transitions. The possibility of spin manipulation by means of an applied bias voltage 
has revealed electron spin to be an information carrier in both classical and quantum 
information technology[15, 16, 17]. A spin-1/2 particle is an appropriate two-level 
quantum system. This forms a quantum bit (qubit) being the basic unit of quantum 
information[18, 19, 20]. 

The full understanding of spin evolution in two-level atoms is necessary to establish 
predictions in AT-level atoms where more complicated patterns are expected as the 
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total spin is «S'[21]. A breakdown between theoretical and experimental expectations 
occurs when the spin dynamics become cumbersome [21, 22, 23]. This is precisely a 
genuine signature of the complexity of spin dynamics in synthesized nanomagnets [24], 
The fundamental Schrodinger equation for such situations involves M coupled linear 
differential equations with nontrivial solutions. This difficulty is overcome by an 
alternative and rigorous approach introduced by Wei and Norman (WN)[25, 26, 27] 
and corresponds to one of the Lie exponential forms. The WN approach involves 
the arrangement of the time-evolution operator (TEO) into a parametrized form that 
depends on the group dimension n. The time-dependent parametric functions for 
SU{2) symmetry group fully determined the TEO and are universal for all irreducible 
representations of the group [28]. 

Group theoretical requirements of SU{n) are useful for quantum optics[29, 30], 
evolution of squeezed states[31], propagation of Gaussian beams[32]. The present paper 
reveals the role of SU{2) group in the physics of AA-energy levels evolving in time with 
the same symmetry and crossing at a single point by linear variation of a controlling 
parameter (time, coordinate, energy, chemical potential, flux etc). The representation of 
the TEO that we consider was first proposed by Kolokolov[33] for functional integration 
of quantum magnets. The application of additional transformations permits to achieve 
the so-called Wei-Norman-Kolokolov (WNK) representation [33, 34]. The case where the 
linear varying parameter is a regular field refers to as the LZ theory [35, 36, 37, 38]. 

The paper is presented as follows: Section 2 presents the TEO, its matrix elements 
(scattering matrix) properties; Section 3 is devoted to application to multi-level LZ 
transitions. Section 4 illustrates the method for the case when an arbitrary spin S is 
subjected to a regular magnetic field and a classical transverse noise with Gaussian 
realizations. Here, the two limits of the noise being fast or slow comparing the 
characteristic time scale of the noise with that of the system are considered. Section 5 
is the summary. 

2. Wei-Norman method for the SU{2) finite dimensional Lie gronp 
2.1. The time-evolution operator 

Let |'0(^)) be the Schrodinger state vector of an arbitrary spin S at time t. The time 
evolution of the coherently driven system is described by the propagator U(t, to) = 
\'ijj(t)){'il:(to)\, (t > to) which satisfies the time-dependent Schrodinger equation {h = 1); 

tj^U{t,to)='H{t)U{t,to), (1) 

and the initial condition U{to,to) = 1 (where 1 is a unit matrix). The Hamiltonian, 
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wW = 5^e,(()S;/, (2) 

i=i 

is Hermitian and for arbitrary t belongs to the finite dimensional Lie group SU (2). Here, 
Qj{t) is a set of linearly independent complex functions of time. The basis for the Lie 
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space is chosen such that for all generators Sj of SU{2) we have [Si, S'j] = c^Si where 
cS are constants structure on SU{2). 

For our case, the fundamental solution U{t,to) = T exp[—i Jl'^'H{T)dT] (with T 
being a time ordering operator) of equation (1) belonging to the group SO {3) is 
transferred to the TEO U{t) selected as the following Lie exponential form[25, 33, 34]: 

U{t) = exp{f{t)S_) exp{h{t)S^)exp{g{t)S+), (3) 

where S± = S^ ±iSy. Substituting the solution in equation (3) into equation (1), the 
functions f{t), h{t) and g{t) are found to satisfy the following system of differential 
equations [34]: 

r if = 

I ih = Q,{t) + Q.{t)f, (4) 

yig = 

assuming the initial values /(to) = h(to) = gifo) = 0. Here, 0±(t) = 0x(t) ±iQy{t); the 
dot on the functions refer to time derivatives. Solutions of equation (4) depend on the 
Riccati equation [39] (hrst equation in (4)) that can be transformed in to a second order 
differential equation [39]. 


2.2. Scattering matrix and properties 


The state vector |'0(t)) of an arbitrary spin S describing 25* + 1 accessible or non 
accessible states in a time-dependent held may be decomposed as a direct sum 

s 

l^(^)) = Cm{t)\S,m). (5) 

m=—S 

Here, Cm{t) are amplitudes for hnding a projection m of S' on a given axis of the basis 
\S, m) {m = —S, —S '-!-!,..., S). This basis can be constructed with the aid of multispinor 
technique[28, 40]. We hnd the matrix elements (scattering matrix): 

= {S,m'\U{t,to)\S,m), (6) 

of the TEO U{t,to). They relate the amplitudes Omit) through the relation: 

s 

(7) 

m'=-S 

We choose the basis to make the states mutually orthogonal {{S,m\S,m') = 6mm') and 
the closure relation J2m=-s IS", m)(S', m| = 1 always satishes. Then, assuming that the 
system evolves from to = —oo, after some elaborated algebra, the matrix elements of 
the TEO dehned in equation (3) are obtained for arbitrary times t as follows: 


(^ - m)!(^ + m')! 


Um,m'iQ) = A / 7^-- 

’ y (S -|- m)!(S — m )l [m' — my. 


( 8 ) 


rrs {S-m')\{S + m)\f'^-^'e^'^ 9 , 

Um',m{Q) = \ - ^7 - jrPm'miQ)^ m < m, 9 

’ y [b + m')!{b — mj! [m — m'j! 
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and satisfy the closure relation Ylim'=-s\^mm'{Q)?‘ — ^ fixed m] Q = {f,h,g). 
Here, m'iQ) ~ 2 Fi{—S + m', S + m' + l,m' — m + 1; —fge^) with 2 Fi{...) being the 
Gauss hypergeometric function[41]. Also, if /, h and g are real, [Umm'ify^ydW — 
U^, m{—g, —h, —f). If one assumes in this case that = {—g,—h,—f), then 
[U^,m'iQ)]^ = Similarly, when m = m', = e"^^w^^{Q). Note 

that the transition amplitudes between the states with extremal spin projections are 
U§g{Q) = since w§g{Q) = 1 and U^g_g{Q) = + e^fgY^. We present 

U^m'iQ) some values of S: 




f/(l/2) _ 

( ge^F 

y je^F + jge^) ^ 

)' 





( g/l 

\f2e^g e^g^ 

\ 




17« = 

A/2eV 

1 + 2e^fg P2g{\ ^ e 

^fg) 

1 

(10) 



e^P 

\/2/(l + e'^/^) e'^Q + fge^yj 



^(3/2) ^ 







/ g3h/2 

y^g3V2^ 

V3e^^Fg2 



\ 



gV2(l _|_ Sfge’'-) 

e^Fg{2 + 3fge^) 

y^g/i/2^2(i _|_ Jg^h 



x/3e3V2/2 

gh/2y^2 + 3fge’^) e 

-'*/2(l + /^g^)(l + 3/^g^) 

y/3e-’''Fg{l + fge^) 

2 


^ g3A/2 j3 

\/2gh/2j2(i _|_ 

\/3g-fe/2j(l _|_ fge^Y 

g-3h/2j']^ _|_ fgPY 

) 



The matrix element in position {i,j) corresponds to the transition amplitude from the 
state with index j to i. 

With the help of the parametrization in equation (3) and the Baker-Campbell- 
Haursdorff expansion for products of exponential operators, we compute 


U-\Q) 


dU{Q) 

dQ 


= A{Q)G{S). 


Here, 




^(Q) = 


V 


0 

0 


-2ge^ 

1 

0 


-9 

9 

1 


2gh \ 


( 11 ) 


( 12 ) 


/ 


and G{S) = [S'-, 5*2, 5'+]^ is a vector whose elements are generators of the Lie algebra 
on SU{2). This is useful to establish some properties of tj{Q) and find the invariant 
integration measure [42] associated with the parametrization (3). Thus, the invariant 
measure dQ = \A{Q)\dfdhdg on the group SU{2) reads 

dQ = e’^dfdhdg. (13) 

This measure is normalized to the ” group volume'' Vgroup = ^^dfdhdg such that 
dQ = e^dfdhdg/Vgroup is a normalized invariant measure. As a consequence, the 
invariant measure dQ permits the normalization of the matrix elements in equations 
(8) and (9) (see also reference [42]): 

^ A A A 


dQ U='(Q) 




2j2 + 1 


(14) 
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By expanding an arbitrary function of Q into the basis of and making use of 

the integral relation (14), it is demonstrated that: 

+1) [ui.JQ)] [uLAQ')] = '5(Q - Q')- (15) 

jmk 

Consider the positive axis. The system of equations (11) contains three coupled 
equations for dlJ/df, dtj/dh and dlJ/dg. After decoupling, the functions tJS_, US+ 
and USz are extracted and sandwiched at the left and right hand sides by the state 
vector \S,m): 


C±Ul^'(Q) = y(STm)(S±m+!)£/**, 


(16) 

( 17 ) 


where 


(18) 

(19) 

( 20 ) 


jC.+ = dg, C,^di,-gdg, £_ = e dj + 2gd^-gdg, 

and Oq — djdQ. For the negative a-axis, we have 

t. = 9 /, £, = -dh + / 9 /, £+ = + 2/4 - f'dj, 

satisfying the following commutation relations 

[£+,£_] = 2£„ [£„£±] = ±£±, 

where £ = (£, £) are isomorphic to the SU{2) group according to (20). On the other 
hand, jC± satisfy similar relations as jC±, except that obeys: 

At/® „.(Q) = -mV® „,(Q). (21) 

One proves that the invariant Casimir operator 

C? = e ^dfg + dhh + dh, (22) 

commutes with all operators in equations (18) and (19) and satishes the eigenvalues 
equation 

i:V®„,(Q) = S(S + l)f/®„,(Q), (23) 

By writing out the /-, h- and ^f-dependence in equations (8) and (9), we obtain the 
differential equations for when m! > m, 


e-^dt 


do + \m' — ml - 
• 9 


+ dh 


dh + [2m' + 1] 


-(.S-m')(^ + m' + l)|M;^^,^,(Q) =0, (24) 


e ^dg df + [m — m'] j 


+ dh 


dh + [2m + 1] 


-(^-m)(^ + m +l)|M;^,„,(g) = 0, (25) 


and 
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when m' < m. The change of variable ^ = —fge^ yields the eqnation for Gauss 
hypergeometric function 2 -Fi(...) respectively related to equations (24) and (25): 


r / N 

r -| 


1 

1 

+ 

1 

to 

+ 

w 


+ (^-m')(^ + m' +1 )|m;^^„,(z) = 0, (26) 


r . 

r -| 

1 

Li- 

R-l , 

O 
to 1 

+ 

1 

1 

+ 

1 

to 

+ 


d 

dz 


+ {S-m){S+ m + l)^wi,^^{z) = 0. (27) 


l-2z 


Relations (22)-(27) are valid for £ as Thus, equation (23) leads to a unique 

equation which holds for all m and m': 

4-S(S + 1)}(7®„,W = 0. (28) 

An alternative form of equation (23) is obtained by changing U^^^,(Q) e~^U^ ,^i(Q). 

This procedure leads us to, 

(29) 


, . d? 


^ ^fg + — dh — S(S + 1) — 0. 


The properties (8)- (29) form a solid background to discuss problems of transitions 
between the 2S' + 1 multiplets of a spin S placed in a time-dependent magnetic 
held. This reminds us of the universal Wigner T)-matrices[43].The algebra of these 
matrices span into the SU(2) group and is frequently used to describe spin dynamics. 
A similar system of three equations like (4) associated with the parametrization of 
the TEO (rotation operator) with time-dependent Euler’s angles involves complicated 
combinations of trigonometric functions with non trivial solutions [44]. The contrast 
with the parametrization (3) becomes evident as the system of equation (4) depends 
on a known equation (Riccati equation) with computable solutions. Some examples are 
shown below for illustration. 


3. Landau-Zener transitions for an arbitrary spin S 

Transitions between energy levels at avoided crossings are known as LZ transitions. We 
write down the TEO for the LZ model generalized to the 25* -|- 1 levels of a spin S. The 
vector 0(f) is identihed by 0Lz(t) = (2A, 0, 2vt) where v is the Zeeman constant sweep 
velocity, A the coupling between split levels. 

For the sake of generalizations of statements and from the 25* -|- 1 dimensional 
representation of the universal enveloping algebra of SU{2), the model is written in the 
basis of the irreducible spin tensors as follows: 

25 L 

L=Q M=-L 


( 30 ) 
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“iflf (^) linearly independent complex fnnctions of time. The matrix elements of 

in the basis \S, m) are calculated with the aid of the Wigner-Eckart[40, 45] theorem and 
are written as follows: 


T; 


( 5 ) 


LM 


2L + 1 


2^ _j_ 


(31) 


where = {Sm', LMWS^m) are Clebsch-Gordan coefficients [40, 45]. The 

linearity of the LZ Hamiltonian and the Wigner-Eckart[40, 45] theorem impose L = 0 
and L = 1 in equation (30). Thus, 


(S) _ 

^0,0 


CTn n — 0, 


= y2A3l(®>, aS5i = -y2AA<®l, aS5 = 2a«3l(®>, 


(32) 


where 


^ ^ S(S + 1X23 + iT 


(33) 


The TEO is fully determined once the Riccati equation is solved. We transform the 
Riccati equation in (4) to a second order differential equation with the transformation 
f(t) = InZ(t) where Z(t) = 3^(^)exp Qz{T)dT^. The function 3^(0 satishes 
the parabolic cylinder (Weber’s) equation[21, 41, 46]: 


de 


3^(0 + 


id 


1 

2 


4 


Here, ^ = y/^te and <5 = A^/n. 
3^(,^) = y{t) as follows [21]: 


3^(0 = 0. (34) 

We finally obtain f{t), h{t) and g{t) in terms of 


A*‘‘'^<*>-A’ 


h{t) = 21ny(f), g{t) = —iA 


dt' 


'to 




.(35) 


The parameters f{t), h{t) and g{t) entirely determine the TEO. They do not contain 
any representation index of the group SU{2). They are universal for all irreducible 
representation of SU{2). This is a generalization of Majorana[14] solutions (solutions 
in terms of parabolic cylinder functions). LZ transitions probabilities depend only on S 
and S. 

Equation (34) has four solutions [41], with an appropriate choice of the solution 
which satishes 3^(—cx)) = 1, 


/(oo) = ye h{oo) = -iid, g{oo) = (36) 

where y = [2sinh(7r5/2)]^/^e’’'^/^ and $ = ^ + argr(l — id) + (5(ln5 — 1) is the Stoke’s 
phase originating from large expansion of parabolic cylinder functions[47, 48]. 

On the other hand, the product of two Gauss hypergeometric functions leads to 
Appell’s function iTi)...] [49, 50]. This permits to arrive the generalized formula: 


piS) 

' m'^m 


{oo) 


{S + my. 

(S'— m')! (m'— m)!(m'+ m)! 


(37) 


for m' > m where 


= F^—S + m', S' + m' + 1, m' — m + 1, m' + m + 1; a^; 5^), (38) 












a = 1 — e~'^^ and b = We achieve our goal. The case m' < m is obtained from 

(37) by swapping m ^ m' everywhere they appear. This compact formula may be used 
to work out the transition probabilities for exceedingly large values of S. 

We present some limiting cases. For non-adiabatic transitions between the states 
with extremal spin projections, V^_lg{oo) = and V^g\_g{oo) = . We achieve 

the LZ formula when S = 1/2. Similarly, for non-adiabatic transitions between the 
state with maximum spin projection and other projections, = 1. The transition 
probabilities are given by 


V 


is) 

S^m 


,S-r. 




oo) = 


{S-m)\{S + m) 


i2S)\. 


(39) 


These probabilities are independent and obey the condition X]m =-5 ~ 

probabilities for transitions between intermediate states can be calculated from (37). 


4. Transverse noise on an arbitrary spin S 

The spin S interacts now with a noisy environment (host lattice; nuclear spins that 
create a held of random amplitude and direction, the Overhauser held associated with a 
classical noise with Gaussian realizations). The transverse components of the magnetic 
held randomly huctuate while its longitudinal component oscillates with a frequency 
u{t) = 2vt: 

0(t) = (40) 

The helds 77* (t) {x, y) are realized as zero-mean Gaussian noise with two-time correlation 
functions given by = 77 jj( 7 |t — t'\) where 

- t'\) = r]^ exp{-'y\t - t'\)6ij. (41) 

Here, rj is the noise strength while 7 is the inverse characteristic correlation time 
of the noise. We discuss noise-induced transitions in the frame of equation (4) and 
construct the relevant Fokker-Planck (FP) equation [51, 52, 53]. Equation (4) however 
contains fast oscillating terms which should be removed. To do that, we introduce slowly 
oscillating variables, 

fit) -t fit) exp (i f 0 ^(r)dr), h{t) h{t) - i [ 0 z(r)dr, g{t) g{t). 

^ Jto ' Jto 

(42) 

The stochastic Langevin variables become Qit) = {f{t),h{t),g{t)). We examine in the 
next sections two limiting cases of noise. 

4 . 1 . Fast noise for arbitrary spin S 

For this case, 7 —?■ 00 , the characteristic time scale of the system is exceedingly large 
compared to that of the noise. With the help of Novikov’s theorem[54], a relevant FP 
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equation is derived for the probability density function (pdf) W(Q,f) as follows(see 
Appendix A for details): 

dtW{Q, t) = -n{t)jCFp{Q)W{Q, t), (43) 


with the initial condition W(Q,to) = S{Q — Qo) where Qo is an initial conhguration; 
Q{t) is the noise spectral density capturing information about the environment (see 
Appendix A). 

^FpiQ) = e~^dfg + 2dfhf - djfp + Off - dh, (44) 


a FP-like operator. From equations (18)-(23), equation (44) may be transformed to the 
following: 

^fp{Q) = ^fp{Q) — i^z + 1)^, (45) 

where 


^FpiQ) — e ^dfg + dhh — dh- (46) 

The operator (£^ + 1) in equation (45) enters the FP operator Cpp as a consequence 
of the non-stochastic character of the magnetic held along the direction of the sweeping 
held. For its random distribution it is found to vanish. Interestingly, W((5, t) in the FP 
equation will be entirely related to It is instructive to note that: 

(4 + (47) 

whereby, Cfp{Q) has the eigenvalues L(L-|-1) —and the eigenfunctions 
After solving equation (43), one obtains 2S'-|-1 solutions. The hnal solution is expressed 
as a linear combination of all particular solutions with 25* -|- 1 constants. The correct 
solution should satisfy the initial condition W{Q,to) = 6{Q — Qq). The normalization 
constants that enters the hnal solution are determined with the aid of property (15). 
The pdf is then obtained as follows: 


W{Q,t)= Y, i‘^L + l)e 


-[L{L+l)-M'^]8{t)/4 


u 


M',Ar(Qo 


-1 r 




LMM’ 


where, 9{t) = 4 Vt{t')dt' is associated with the phase accumulated in an interval (to, t]. 
For a sweeping varying from to = —00 to t = 00 , one hnds that 6{oo) = 47r77.(0)/u (see 
Reference [55] for more details). 

We compute the probability for an arbitrary spin S initially in a conhguration 
Qo at time to to end up at a diherent conhguration Q at time t by expanding the 
density matrix into the basis of the irreducible spin tensors (analogous to spherical 
harmonics) [56]: 

2S L 

'5W = E E FhlXiuS- («) 

L=0 M=-L 

The tensors operators form an orthogonal basis in the (2S'-|-1) x (2S'-|-1) space of the 
density matrix which corresponds to SO{3) rather than SU{2) group. The expansion 
(49) is valid for any Hermitian matrix with a trace equals to 1. The brackets (...) express 
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Tr(...p), the coefficients of expansion (49) are Bloch tensors. They evolve 

independently for all given values of M according to an equation which is analogous to 
the optical Bloch equation for a Bloch vector set on a unit sphere. (T|'^)(f) can also 
be viewed as an average of over environmental [env) degrees of freedom when 

P —t peni; (reduced density matrix). For our case, (...) accounts only for spin fluctuations, 

(...) -)■ {...)phase with 

{...)phase{t) = j dQ...W{Q,t). (50) 

The distribution W(Q,f) is spatially uniform (no spatial decoherence). We hnd 
the Bloch tensors for S in the initial conhguration Q = Qq. The procedure in equation 
(50) yields 

(fl2)(t) = 


The diagonal components of the Bloch tensors (51) are non vanishing only for zero 
projection [M = 0) of S' onto the direction of the sweeping field. The non diagonal 
components are different from zero if there is a coherence in the initial state. For a 
different projection (M ^ 0), the Bloch tensors vanish as a result of averaging over 

/' / o\ 

the random phases in the initial state. The diagonal components (T^q/(^) Bloch 

tensors involve Sz, S'^S' 2 ,...They are combinations of diagonal components pjj of the 
density matrix (transition probabilities) and represent population differences between 
levels. For instance, = ^/2{pll{t) — p 22 {t)) is the population difference for 

two-level systems. A key observation related to (51) is that the diagonal components 
(T^o/(^) of fho Bloch tensors decrease monotonously with time due to the noise. 
This remark shows that the population differences in average can only decrease after 
transitions. They may increase if there is a coherence in the initial state as non diagonal 
components of the Bloch tensors are not zero. 

Considering equation (15) and (49) then, the transition probability is obtained as 
follows: 



1 

2F + 1 


2S 




2L + 1 
2F + 1 


rp(S) 

-‘-LO 


^S,m' L{L+l)e/A 


.(52) 


Similar results were obtained in Reference [58]. When S' = 1/2 and S' = 1 we reproduce 
the results in references [55, 57, 58]. For illustration. Tablet shows the transition 
probabilities between Zeeman multiplets when S = 3/2. The remaining part of the 
table being completed using symmetry between levels or identical procedure. Equation 
(52) is obtained with the help of [56, 58]. 


4-2. Slow noise for arbitrary spin S 

For the case of slow noise, 7 —)■ 0, the correlation time is exceedingly greater than 
the characteristic time of the LZ process. A FP equation is not easily achieved. We 
readily solve the Langevin equations in phase space {f,h,g). We consider the long¬ 
time asymptotic transition probabilities when t ^ 00 and perform Gaussian ensemble 
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State-to-state Final occupation ioi t = oo 


3\ V |3 3\ 1 

2/ 12’ 2/ 4 

1 + ^ ^-36/2 ^ |g-0/2 

3\ ^ |3 1\ 1 

2/ ^ 1 2’ 2/ 4 

1 _ 3g-30 _ g-36(/2 

3\ ^ 13 _1\ 1 

2/ 1 2’ 2/ 4 

1 p Iq-36 _ g-36»/2 _ |g-0/2 

3\ |3 3\ 1 

2/ 1 2’ 2/ 4 

1 _ l(.-3d p g-3e/2 _ |g-e/2 

-1\ ^ |3 1 

2/ ^ 1 2’ 2/ 4 

1 p |e-3e _ q-39/2 _ |g-e/2 

_1\ ^ |3 1\ 1 

2/ 1 2’ 2/ 4 

1 _ |e-3^ p £-36/2 _ 

_1\ ^ |3 _1\ 1 

2/^12’ 2/ 4 

'l P |e-3^ p e-3G2 P ie-^/2 

_1\ ^ |3 _3\ 1 

2/ 12’ 2/ 4 

1 _ |e-3® _ ^-39/2 |g-6/2 


Table 1. Landau-Zener transition probabilities for some of the four multiplets 
of the spin S = 3/2. (see also references [57, 58]) 


average {...)noise (given in reference [55]) over the distribution Q of the noise[55] i.e.: 





(53) 


For a single component transverse noise (X-noise, (3 = X), Q = x while for two- 
component (XF-noise, (3 = XY), Q = \/aF+^ and the average {...)noise is transferred 
to a double Gaussian distribution (see reference [55]). For transitions between the state 
with maximum spin projection {m' = S) and that with other projections (m < —S), 
the procedure (53) accounting for (39) leads to 


S—m 


^ = E u.(s + m)\ 


1)^(2^)! 




m 


L)\V 


l + ^(^ + m + L) E(54) 


Here, k takes two values 1 and 2 for single- and two-component transverse noise 
respectively. Some results illustrating this procedure are presented below for S = 1/2, 
^ = 1 and ^ = 3/2: 


p(l/2) 


P Q 
q p 


p(i) 

sn 


p(3/2) 

sn 


^ p^ 2pq q^ 
2pq {2p — 1 )^ 2pq 
^ q^ 2pq p^ 

( p3 3p^q 

Sp'^q (3p —2)^p 
3pg^ (3p — l)^g 
y q^ 3pq^ 


\ 


/ 

3pq^ 

{3p - l)2g 
{3p-2fp 
3p^q 


\ 

3pq^ 

3p^q 

p3 y 


(55) 

(56) 


( 57 ) 
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Here, q = 1 —p with = [1 + ]~t. It is instructive to note that when using q and 

p we hrst collect all powers of p. For instance pq = p — p^ oi pq^ = P — + p^. 

Slow noise transforms the Gaussian properties of LZ transitions probabilities to 
Lorentzian[55]. 

5. Summary 

We have developed an approach to problems of transitions in AA-level systems induced 
by linearly varying protocols (time, coordinate, energy, chemical potential, flux etc) 
and(or) fluctuations of Overhauser held associated with hyperhne interactions. The 
approach is applied to cases of transitions induced by random transverse (or regular) 
magnetic helds with focus on level-crossing systems with SU (2) symmetry. The universal 
model Hamiltonian that considers all irreducible representations of an arbitrary spin S 
is achieved by expressing the Landau-Zener Hamiltonian through spin tensors which 
envelop the universal su{2) algebra. 

The fundamental solution to the time-dependent Schrodinger equation is expanded 
as a product of three exponential operators (with three time-dependent functions /(t), 
h{t) and g{t)). We indicate the condition under which the TEO is fully determined, 
namely once the function f{t) (solution of a Riccati equation) is obtained. For general 
SU{2) Landau-Zener problem, this equation is analytically solved, /(t), h{t) and g{t) 
are expressed through Weber’s parabolic cylinder functions. The huger prints of SU (2) 
irreducible representation indices disappering. This is observed to be a generalization 
of Majorana solutions[14] for two-level crossing. A compact formula expressing LZ 
transitions between diabatic states is derived through the Appell function F4(...)[49, 50]. 
The expressions of the SU{2) ’’scattering matrices” and their useful properties are 
presented. 

Considering noise-induced Landau-Zener transitions, the set of equations for f{t), 
h{t) and g{t) are observed to be described by nonlinear Langevin equations. In the 
limit of fast noise, a relevant Fokker-Planck equation for the probability density function 
W((5,f) is derived. The properties of the scattering matrix elements U^^,{t) help to 
solve this equation exactly in phase space (/, h, g) where spin dynamics are described. 
The originality of the approach presented in this paper consists in its generalization. 
Rather than solving the von-Neuman equation for the density matrix we have expanded 
the latter in the basis of irreducible spin tensors. The coefficients of this expansion 
(Bloch tensors) were obtained by analytically solving a Fokker-Planck equation. In 
the expression for probabilities of diabatic transitions between Zeeman multiplets we 
introduce functions with clear physical meaning such as Bloch tensors. We express 
our results in terms of well-know mathematical functions such as the Appell function 
F4(...)[49, 50]. It is observed that all transition probabilities are exclusively tailored by S 
and(or) /v{rf‘/v). The results for fast noise are observed to be in excellent agreement 
with those in reference [57, 58]. In the opposite limit of slow noise, the nonlinear Langevin 
equations are solved and a Bloch average procedure [55] is adopted. 
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The present study helps in explaining the mechanism of spin-flip transitions and 
spin-propagations [4] in single- and multi-molecule magnets, synthesis and magnetic 
measurements of cubic nanomagnets [5, 6, 7, 8], quantum computation [9, 10] and 
Landau-Zener transitions [11, 12, 13, 14], The study might also be extended to different 
types of magnetic helds (not necessary linearly varying) and lead to solving a nonlinear 
Ricatti equation [39]. 

Appendix A. Fokker-Planck equation 

We write a Fokker-Planck (FP) equation for an arbitrary Gaussian process in the limit 
of fast noise. We use Einstein’s notation. We consider Langevin[51] equations derived 
from TEOs for the group SU{2). Three generators are involved and result in a three- 
dimensional Langevin equation of the form[51] 

Qtiit) = t) + B^kiQi (A.i) 

yU = 1,2, 3. For the case when the particle is placed in a magnetic held random in the 
transverse directions and regular in the longitudinal direction with a frequency Qz{t), 
the equation (A.l) contains fast oscillating terms that are eliminated by dehning an 
appropriate gauge transformation. The gauge in such cases, consists of tracing out 
Qz(t) in A^{Q,t) (see equation (42)). After doing this, we further consider only slow 
variables Q(t) and omit the sign tilde on Q{t). The noise helds Jx,y{t) are transformed 
as 

Jx{t) = Jx{t)cos'd{t) + Jy{t) sill'd(t), (A.2) 

Jy{t) = Jy{t) cos'd{t) - Jx;{t) sint9(t), (A.3) 

where 

■dit) = j dTOzir), (A.4) 

Jto 

is the phase accumulated by the system during a sweep. As a consequence, noise helds 
remains centered {Jx,y{t)) = 0 but multiplicative: the transverse components of the new 
noise helds are statistically dependent. Equations (A.2) and (A. 3) show that the gauge 
which eliminates fast oscillating terms rotates the coordinates frame of the transverse 
plane of an angle 'd{t). 

We construct a FP equation for the probability distribution function[51] 

W(Q,t) = (5(Q-g(t))), (A.5) 

which satishes the initial condition W(g,fo) = — Qo). The brackets (...) denote an 

averaging over all possible realizations of the noise. If one writes the equation for the 

density S(Q — Q(t)) of the system in the conhguration Q(t) at time t and average the 

resulting equation over the noise realizations[51]: 

dW d d ~ 

^ ^ - Qit)))- (A.6) 
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The remaining average in this equation is simplihed with the help of Novikov’s 
theorem [54] and after some simplihcations yields 

{MmQ - Qm = - Qm. ( at ) 


dQv 


where 


Jku{t)= dtiOkmillt-ti\) 

Jto 


SQu{t) 


(A.8) 

SJmitl) 

and Gkmillt — ti|) are two-time correlation functions for the noise helds J'l, that can be 
derived from (A.2) and (A.3) considering (41). Due to the random character of noise 
in the transverse direction, Jz{t) = 0 and J^3u(t)=0. The relation (A.8) is valid for 
arbitrary zero-mean noise held with Gaussian realizations. 

When the noise is fast, its characteristic correlation time I/7 is small compared to 
the transition time of the system. We consider the noise contribution at the vicinity of 
ti ~ t. The Taylor series expansion of 5Qy{t) j5 Jmiti) is considered. 


5Qy{t) 6 Q^{t) 


dJmiti) dJmiti) 


d SQ,y{t) 


tl=t 


dh 


ti=t 


T 0[r^]. 


(A.9) 


Here, r = f is a characteristic time interval in which the coherence of the system is 
preserved. On the other hand r is small enough as 7 —)■ cxo. This permits to consider 
only the hrst term in the expansion (A.9). We readily integrate the Langevin equation 
to hnd this term and we obtain 


Q.{t)= / dt'A,{Q\t')+ / dt'B,k{Q\t')Jk{t'), 


(A.IO) 


'^0 


/£o 


since Qu{to) = 0. The functional derivative of equation (A.IO) yields: 


dQAt) ^ ^ 


dJmitl 






Hi 


dJm{t 


1 


where 


dQy{t) dQy{t) 


dJmitl) dj'rnitl) 


t=tl 




(A.ll) 


(A.12) 


The integral in (A.8) can be evaluated. To do this, we hrst calculate all the two- 
time correlation functions Qkm{---) in equation (A.8). From equation (41), we write the 
function 7^(7|f — t'\) = ?7^exp(—7|f — t'\) that permits the calculation of the following 

poo 

^ Bi,i{Q{t),t) / dr7^(7|r|) cos[?9(f) -'d(fr)]-h 

Jo 

poo 

Bu2{Q{t),t) dTTl{'y\T\)sm['i}{t) -'d{t + t)], (A.13) 


and 


^-B^i{Q{t),t) / dTTl{'j\T\)sm[d{t) - 'd{t + t)] + 


B^ 2 {Q{t),t) / dr7^(7|r|) cos['d(f) -'d(f-f r)], (A.14) 
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As r is small, the phase 'd{t + r) is expanded as a series of r up to the linear term 
+ r) 'd(t) + + C>(r^). Therefore, we deduce 'd(t) — + r) —uj{t)T where 

the frequency u(t) reads 

= ^ = (A,i5) 

From equations (A. 13) and (A. 14) we have 

Jku{t) ^ (A.16) 

where 

poo 

f2(t) = / dr7^(7|r|) cos[ci;(t)r]. (A.17) 

Jo 

Finally, a FP equation associated with the noise is written as: 

<9W d d‘^ 

where the drift and diffusion coefficients are respectively given by 

D,(Q,t) = A,(Q,t) + (A.19) 

and 

Dfj,u{Q,t) = B^k{Q-,t)Jkv{t). (A.20) 

For our case, the procedure for elimination of fast oscillating terms in equation (4) are 
indicated in equation (42). Q^{t) = [f{t),h{t),g{t)], = Jx,yif) and JziJ) = 0. 

The column vector A^{Q, t) vanishes while the elements B^k{Q-, t) form a (3 x 3) matrix. 
Qz{t) = a;(t) = 2vt. From equation (A.18), after a long calculation, we derive the FP 
equation in (44). 
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